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ABSTRACT 

The population of compact massive galaxies observed at 2: > 1 are hypothesised, both 
observationally and in simulations, to be merger remnants of gas-rich disc galaxies. 
To probe such a scenario we analyse a sample of 12 gas-rich and active star forming 
sub-mm galaxies (SMGs) at 1.8 < z < 3. We present a structural and size measure- 
ment analysis for all of these objects using very deep ACS and NICMOS imaging in 
the GOODS-North field. Our analysis reveals a heterogeneous mix of morphologies 
and sizes. We find that four galaxies (33%±17%) show clear signs of mergers or 
interactions, which we classify as early-stage mergers. The remaining galaxies are 
divided into two categories: five of them (42%±18%) are diffuse and regular disc-like 
objects, while three (25%±14%) are very compact, spheroidal systems. We argue that 
these three categories can be accommodated into an evolutionary sequence, showing 
the transformation from isolated, gas-rich discs with typical sizes of 2-3 kpc, into 
compact 1 kpc) galaxies through violent major merger events, compatible with the 
scenario depicted by theoretical models. Our findings that some SMGs are already 
dense and compact provides strong support to the idea that SMGs are the precursors 
of the compact, massive galaxies found at slightly lower redshift. 

Key words: galaxies: evolution — galaxies: high-redshift — submillimetre — galax- 
ies: starburst — galaxies: active 



growth in galaxy sizes while permitting a mild evolution 
in velocity dispersion as ob served (e.g. iCenarro Truiillol 
I2OO9I : ICapDellari et aLllioOol ). 

An additional problem is understanding how these ob- 
jects first formed. An em erging pictu re (e.g.. [ Hopkins et al.l 
I2OO7I : ICimatti et aDIioOS : Hopkins et al]l2009d ) for the for- 
mation of these compact systems predicts that massive, gas- 
rich galaxies at very high redshift become unstable follow- 
ing a major merger event, triggering a sho rt-lived starburst 
within O.lGyr. Theoretical models (jKhochfar Silkl 
I2OO6I : Ih opkins et al. I I2OO7I I have shown that the size of the 
remnant strongly depends upon the degree of dissipation in- 
volved, being very small in the case of strongly dissipative 
mergers. Since at high redshif t (z > 2) gala xies are more 
gas-rich than they are today (|Erb et al.ll2006l l. the degree of 
dissipation is expected to be high, and the resulting remnant 
extremely compact, with sizes < 1 kpc. 

Given the great amount of gas involved in these star 
formation processes, we expect the progenitors of mas- 
sive compact galaxies to be undergoing a high amount 
of star formation, and hence should be detectable in 



1 INTRODUCTION 



Observations have shown that massive spheroids at high 
redshift are rem arkably smal le r than their local counter- 
parts (e .g. Da ddi et al.1 l2005l: iTruiillo et al.1 E)0 6a. 200l 
Longhetti et al.1 l2007l: iBuitrago et al.1 12OO8I : ICimatti etld] 



20081 : Ivan Dokkum et al. " l2008l V This population of galax- 



ies must evolve into local massive ellipticals based on 
their stellar masses, however it is not yet clear what 
the primary mechanism is which increases their sizes 
by a factor of four to match the local galaxy popu- 
lation. Although various effects (i.e. dry major merg- 
ers, observational biase s, selection effects, see for instance 
IValentinuzzi et al."201Q') may explain the observed size evo- 
lution of these systems, a physically motivated favoured 
mechanism is the growth in size by later minor mergers 
with less dense galaxies | Bournaud, Jog, Combes '2007"; 
iNaab. Johansson, k Ostr iker 2009; Hopkins et al. 2009a). 
This scenario has the advantage in that it facilities the 
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the subm-mm (iNaravanan et alJ I2OIOI ). To tests this hy- 
pothesis we examine in this paper a sample of SMGs 
(sub-mm galaxies), which have been imaged with HST 
in deep exposures, to probe their structural properties. 
SMGs are among the most luminous (L 1O"'^^L0), 
and rapidly star forming ( SFR ^ 10^Mc;)/igyr ) g;alaxies in 



the hieh-redshift universe (IHuehes et al. 19981: lEales et al. 


2000': 'Smail et al. 2002: CoDDin et al.1 


2006: iTacconi et al. 


2006. 2008: Menendez-Delmestre et al. 


2OO9II and are be- 



lieved to be the prec ursors of local early-type galaxies 
(jSwinbank et al.l I2OO6I ). As their local counterparts, the 
ultra- luminous infrared galaxies, the origin for their high 
fluxes is thought to be a strong starburst and/or AGN activ- 
ity, l ik ely triggered by a ma j or merger (ISanders &: Mirabell 
Il996l : iMurphv et all Il996l : IClements et al.1 Il996l l. It is 
also plausible that many sub-mm galaxies are simi- 
lar to the hy per-luminous in frared galaxies (HLIRGS, 
LiR > 10^^ LrT^. lRowan-Robinsonll200Q ). containing a simi- 
lar a mount of molecular ga s, and an extreme star formation 
rate (jParrah et al.ll2002al lbh. 

Furthermore, the similarity between the stellar mass 
surface densities of SMGs and the comp act massive galaxies 
at lower redshifts (|Tacconi et ahlbood . I2OO8I : ICimatti et al.1 
[2008. ^ makes sub-mm galaxies natural candidates for being 
the precursors of these compact galaxies. This paper tests 
this hypothesis and investigates whether the proposed theo- 
retical scenario for the formation of massive compact galax- 
ies is compatible with present deep imaging and star forma- 
tion analysis of these sub-mm systems. To explore this, we 
probe the sizes of the SMGs in various phases, which we con- 
struct based on morphology and structure, and test whether 
the star formation rates of these galaxies are in agreement 
with theoretical expectations. Second, we probe whether the 
different observed morphologies of the SMGs can be fitted 
into the evolutionary sequence proposed by current massive 
galaxies formation models. 

The paper is structured as follows. In Sect. [2] we describe 
our sample. In Sect. [3] we explain the method adopted for our 
analysis. In Sect. [J] we present our results and discuss their 
implication in Sect.O Throughout, we assume the following 
cosmology: Hq—70 km s~^ Mpc~^, Vim — 0.3 and VLk — 0.7 
and we use AB magnitudes. 



2 DATA 

Our target selection is based on the following criteria, needed 
to characterize the nature of our high redshift sub-mm galax- 
ies: (i) very deep images (particularly with ACS on the 
Hubble Space Telescope) to detect signatures of interac- 
tions; (ii) systems within and similar to the redshift range 
where the massive compact galaxies have been detected (i.e. 
2 < 2; < 3), ideally based on spectroscopic redshifts; (iii) 
our SMGs should be massive (M > 10^^ M©) to ensure that 
these objects are indeed the progenitors of the lower redshift 
massive galaxy population. 

For the above reasons the number of SMGs that 
we can explore is limit ed. We adopt as our parent 
sampl e the SMGs from iMichalowski, Hiorth, WatsonI 
(I2OO9II. This paper c ontains spectroscopic redshifts (from 
IChapman et al.ll2005l l. multi-wavelength photometric data- 
points collected from the literature, including Spitzer ob- 
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Figure 1. NICMOS images of our sample of SMGs with 
ACS contours overlapped. Contour levels of 24,23,22,21 and 20 
mag/arcsec^ are shown with increasing thickness. For each galaxy 
we indicate redshift, logarithmic stellar mass (in units of solar 
masses), effective radius (in kpc) and Sersic index. The white line 
in the lower-right corner is 1 arcsec in size. 



servations. This allows for the determination of accu- 
rate SFRs and ste llar mas ses f or 76 massive SMGs (see 
iMichalowski, Hiorth. WatsonI (|2009l ) for these values and 
the details in how they were calculated) . The resulting range 
of stellar mass is quite small: 10^^ — lO^^M©,. For these rea- 
sons this sample of sub-mm galaxies are the progenitors of 
a fraction of local massive ellipticals, as well as the compact 
galaxies at high redshift. Since we are interested in the high- 
redshift population, we restrict our analysis to the redshift 
range: 1.8 < z < 3. The original stellar masses from this cat- 
alogue are converted fro m a Salpeter In itial Mass Function 
(IMF) to a Kroupa IMF (|KrouDall200lh . 

Our new analysis is based on deep archival HST ACS 
and NICMOS images of these sources, where we find im- 
ages for 12 sources in the GOODS North field. We also 

investigated sources in other fields (e.g., Webb et al 20031 : 

IClements et al.|[2004l l, but these have too shallow observa- 
tions to be useful for our purpose^]. 

Our final collection of images is as follows: we have ACS 
images in the F850LP filter (z-band) for all 12 objects, and 
NICM OS NIC-3 data fr om the GOODS NICMOS Survey 
fGNS: iBluck et al.ll2009l : Conselice et al. 2010 in prepara- 
tion) in the F160W filter (H-band) for 6 objects. The z-band 
data, at 5 orbits depth, reaching a magnitude limit of z=27 
mag (15cr), have been drizzled to a scale of 0.03" with a 
Point Spread Function (PSF) Fuh Width at Half-Maximum 



^ These observations typically have integration times < 7000s 
compared to the > 27000s for the GOODS-N imaging 
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of 0.1". The NICMOS images, at 3 orbits depth (hmiting 
magnitude of H=26.8 mag, 5cr ), were combined to produce 
images with a pixel scale of 0.1" and PSF size of 0.3". Note 
that at z ~ 2.5 the ACS data trace the near-UV rest-frame, 
whereas the NICMOS imaging shows the B-band rest-frame. 
The NICMOS and ACS images of the final sample are shown 
in Figure [1] and [21 respectively. 

All the objects in our sample show hints of AGN ac- 
tivity either from optical spectral features, or from X-ray 
and radio detections. However, it is unlikely that these AGN 
dominate the bolometric emissivity in such objects. The rea- 
son for this is that all the SEDs are well-fit by pure star 
forming models, without the need of an AGN contribution 
(|Michalowski. Hiorth. k Watsonllioool ) . The mid-IR spectra 
of these sub-mm galaxies further confirm that they are star- 
burst domin ated without a bolometrically significan t AGN 
( Pope et"aD 2008). Moreover, as shown by Alexander et al.l 
( 20051 ), the X-ray to far-IR luminosity ratio is much lower 
than that found in QSOs, indicating that the contribution of 
the AGN to th e total luminosity can not be higher than 10%. 
More recently. Laird et al. I (|201 0') found that the X-ray emis- 
sion in SMG is largely due to star formation activity and, 
even in the cases where the presence of AGN is confirmed, 
it is not the dominant contributor to the bolometric lumi- 
nosity, except in rare cases. Therefore we expect that the 
properties of our sub-mm selected galaxies are not biased 
by AGN activity. 



3 ANALYSIS 

Morphological parameters and sizes are measured on our 
sample through t he two-dimensional fitting code GALFIT 
(IPeng e t al.ll2002l ). We model the light distribution of the 
sources with a Sersic profile, deriving the Sersic index n, 
the axis ratio b/a and the semi-major effective radius Ve 
in arcsecs. We scale the radius to the physical scale in kpc 
relative to the redshift of the source and we circularise it 
through Re — re^/hja. The robustness of GALFIT in re- 
covering sizes and structural parameters usi ng HST data 
was assessed in several previous papers (e.g. iTruiillo et al.l 
l2006bl . I2OO7I : ICimatti et al]l2008l l. through the use of simu- 
lated galaxies, and is found to be robust at the resolution 
and depth of our HST imaging. 

To check the reliability of our measurements against 
changes in the PSF shape along the images we use up 
to five different natural stars, as found in our fields, as 
PSFs. We measure the structural parameters for each ob- 
jects five times, taking as our final measurement the bi- 
weight estimator, and its confidence interval as the uncer- 
tainty in the measurement. Most of the sources have an 
uncertainty of < 20% in the measured size and < 30% in 
the Sersic indexT'Only two cases (SMMJ123606. 72+621550. 7 
and SMMJ123632.61+620800.1) contain large errors based 
on the ACS imaging analysis, as their effective radii are close 
to the size of the PSF. 

In the NICMOS images, the most uncer- 
tain measures are found for the compact source 
SMMJ123632.61+620800.1, and for the merging sys- 
tem SMMJ123616. 15+621513.7, where the NICMOS 
resolution does not resolve the three components seen in 
the ACS observations. In most of the cases, contaminating 



neighbours are present in the image, and we account for 
these by fitting the surface brightness profile of these 
neigh bours tog e ther with the target object. As shown by 
Haus sler et al T (|2007l l the simultaneous fit allows one to 
recover structural parameters more reliably than deriving 
them using masks. Likewise, in the case of mergers, the 
interacting systems are fit at the same time. The results 
of our analysis are presented in Table [1] where size, Sersic 
index and axis ratio are shown for the ACS and NICMOS 
data. 

By comparing sizes and Sersic indices measured in dif- 
ferent filters, we find that the two estimates for the isolated 
systems agree to within 3cr. The reason for the slight differ- 
ences in these fits may be ascribed to the morphological k- 
correction, where different aspects of the stellar populations 
for these galaxies are observed at different wavelengths. This 
is particularly relevant within actively star forming galax- 
ies or those with dust, such as the systems we examine in 
this paper. Note that at the median redshift of our sample 
2.3), the NICMOS filter is probing the rest-frame opti- 
cal, while the ACS imaging matches the rest-frame NUV. 
Hence the bluer ACS band is more sensitive to the clumpy 
distribution of the star forming regions. Particularly in the 
case of 'disc'-like galaxies, the Sersic index measured in the 
NUV rest-frame is smaller compared to the rest-frame op- 
tical wavebands (see also Rawat. Wadadekar . De Mellol 
l2009h . However, in most of the cases the NUV and optical 
rest-frame morphologies agree in their ability to discrimi- 
nate between disc-like (n < 2.5) and early-type (n > 2.5) 
galaxies. 

On the other hand, for the merging systems the mea- 
surements in the NICMOS and ACS bands show signifi- 
cant differences. One reason for this is the lack of reso- 
lution in the NICMOS band, which does not always al- 
low us to resolve the merging systems, such as the case 
of SMMJ123616. 15+621513.7 where three galaxies are seen 
in the z-band image, but only two in the H-band im- 
age. The system SMMJ123622.65+621629.7 is a rapidly 
ongoing merger which is seen in the rest-frame UV as a 
compact, elongated galaxy with prominent tidal features, 
while in the NICMOS band the tides are not resolved, and 
the effective radius is much larger. In the case of system 
SMMJ123635. 59+621424.1 the strong disagreement is due 
to the fact that in the z-band we observe a multicompo- 
nent system, made up of a compact bulge with a double 
nucleus, plus faint spiral arms on larger scales. Since the 
spiral arms are too faint to be fit with a single Sersic model 
we have measured the compactness of the bulge. Rather, in 
the NICMOS H-band the multicomponent structure is not 
resolved, and fits to the system reveal a much more extended 
profile. Hence, merging systems can appear very different at 
different wavelengths when using images at different reso- 
lutions, and morphological k-corrections can strongly affect 
our measurements. 



4 RESULTS 

In Figure [2] we show the ACS images of the sample, clas- 
sified according to their morphological properties. We find 
that SMGs display a heterogeneous mix of morphologies. 
We divide the sample into three main categories: disc-like 
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Figure 2. ACS images of our sample of 12 SMGs. The galaxies are ordered according to the evolutionary sequence described in the 
Sec. 4. The first five, from left to right and top to bottom, are the disc-like objects, the next four are the ongoing mergers, and the final 
three are our classified compact systems. Each panel shows a box-size of 40x40 kpc. The white line in the lower-right corner shows the 
1 arcsec scale-length. In the cases where the target object is not obvious it is indicated by two white bars. As in Figure [T] stellar masses 
are in units of solar masses and effective radii in units of kpc. Specific star formation rates (SSFR) are in Gyr"-*^. 



objects (5/12), ongoing mergers (4/12) and compact galax- 
ies (3/12). The first class includes the majority of the sources 
(42%±18%), which show quite regular and diffuse morpholo- 
gies, light-profiles characteristic of late-type galaxies (i.e. 
n < 2.5), and are very faint. Although we can not rule out 
that some of them are interacting systems, where the com- 
panion is not detected due to its faintness at NUV wave- 
lengths 0, our hypothesis is that these objects are disc-like 
galaxies with ongoing SFR. It is important to note that for 
galaxies in this category, where we have NICMOS imaging, 
the disc-like nature is confirmed through the rest-frame op- 
tical images. 

The second category, ongoing mergers (33%±17%), are 
systems where two or more components are clearly visi- 
ble in the ACS band. We include in this class the source 
SMMJ123635. 59+621424.1, since a double nucleus is visible, 
and it is likely to be in the final stages of a merger. Most of 
these interacting systems display disc-like structures with a 
large range in sizes. This is also confirmed in NICMOS when 
data is available. The final class includes three (25%=bl4%) 
very compact, isolated sources. All of these objects show 
a concentrated light-profile, with Sersic indices indicative 



2 Indeed two of the disc-like objects, SMMJ123711. 98+621325. 7 
and SMMJ123707.21+621408.1, present multi-component radio 
counterparts lying at the same redshift (|Swinbank et al.l l2004l : 
IChapman et al.ll2005l ). Therefore, although they appear isolated 
in the optical they a re likely early stage mergers (,Pope et al. 200^ : 
iTacconi et"ani2008l ). 



of an early-type morphology {n > 2.5) and Re < I kpc. 
Note that the extremely small sizes measured in the ACS 
data (tracing the NUV) indicates that the star formation is 
extremely concentrated in the centre of these objects. Al- 
though the AGN contribution might play a role in shrinking 
the size measured in the rest-frame UV, we find that when 
NICMOS data are available the near-IR sizes are slightly 
larger but still very small compared to galaxies with similar 
masses in the nearby universe. 

In Figure \3\ we plot the stellar mass-size relation for 
our sample galaxies, divided into the three morphological 
classes. O verplotted is the mass-size relation for disc galaxies 
from B uitrago et al.l ([2008i ) at z=2.3 (the median redshift of 
our sample) : 

Re{M,,z = 2.3) = a{l + z)-^Re{M,,z = 0.1) , (1) 

with a = 1.1 and (3 = 0.8 for NICMOS-derived sizes (upper 
envelope of the shaded region) and a = 1.1 and /3 = 1.0 for 
the ACS-derived sizes (lower e nvelope). The loca l relation, 
Re{M^,z = 0.1), is taken from ' Shen et al.1 ([2003') based on 
SDSS data. Remarkably, the sizes of the disc-like galaxies 
in our sample match the sizes found for the general mas- 
sive g;alaxy populatio n at this redshift. As the relation of 
iBuitrago et al.l (|2008l ) does not change significantly for ACS 
and NICMOS data, this reinforces the idea that the ACS 
sizes of our disc-like objects (around 2-3 kpc) are represen- 
tative of the rest-frame optical sizes of these galaxies. 

We note that some of the galax ies in our sample were 
already studied in previous works (IChapman et al.l 2003; 



already studied m previous works (luhapman et al.l zUUd : 
IConselice. Chapman, Windhorst l2003l : ISmail et al.|[20o1 : 
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Figure 3. Stellar mass-size relation for our sample. Red squares 
indicate the disc-like galaxies, green circles are for the 'mergers', 
and blue asterisks represent the three compact objects. Large 
and small symbols are for NICMOS and ACS data respectively. 
Overplotted is the stellar-mass relation for disc galaxies evolved 
to z=2.3 (yellow shaded region: upper envelope is the NICMOS- 
derived relatio n, lower envelope is the ACS-derived one), following 
lBuitrap;o et a l. (2008), and the local relation (black dashed line) 
from lShen^ al.. (,2003,). 
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Figure 4. Specific star formation rate for our 12 SMGs (points, 
solid line). Different symbols refer to different morphological 
classes as in Figure [3] (red squares: disc-like galaxies, green cir- 
cles: mergers, blue asterisks: compact galaxies ). The open sym- 
bols with errorbars indicate the average SSFR for a given class. 
Average and standard deviation have been computed with a boot- 
strap resampling. The x-axis refers to the position of the galaxy 
in the proposed evolutionary sequence shown in Figure [2] The 
dashed line shows the predictions of Naravanan et al. (2010) from 
a merger model for SMGs, with the time-scale shown in the upper 
axis in units of Gyrs. 



lAlmaini et aD l2005l : IPodg et al.1 l2005l l using both ground- 
based and HST data. They found a higher merger fraction 
(^ 50 — 60%) with respect to ours 30%), partly due to 
the method adopted for the classification of mergers. While 
the above authors use schemes based on the asymmetry pa- 
rameter, in our analysis we only considered ongoing mergers 
those systems where two or more distinct components are 
resolved in the ACS image. As mentioned above, it is also 
possible that dust is hiding the interacting structure, and we 
are missing a fraction of actual mergers, hence our measured 
merger fraction has to be considered a lower limit. 



5 DISCUSSION 

Motivated by hi^h - resolution hydrodynamical simulations 
(iDekel Co3l2006l; ICox et al.|[200 8b: Hopkins et al. 2 0093 : 
iNaravanan et al.ll2010l l. we argue that the three morpholog- 
ical classes outlined above could represent an evolutionary 
sequence, where the disc-like class represents a pre-merger 
phase followed by the major merger event, while the compact 
sources can be interpreted as the end-stages of the merger, 
caught during or just after the coalescence. Numerical simu- 
lations naturally explain how the tidal forces involved in the 
interactions remove angular momentum from the systems, 
allowing the gas to fuel towards the centre. At this stage, 
the gas is compressed into a very small volu me, leading to 
surface densities of the order ~ 10^ Mopc"^ ([Hopkins et al.1 
|2009b|), close to that of molecular clouds. According to the 
Kennicutt law, the SFR in such conditions is extremely en- 
hanced. Since the dynamical time-scale that drives the col- 
lapse is similar to the star formatio n rate, the system rapidly 
exhausts its ^as whi le it contracts (jMihos k, Hernqu'istll 19961 : 
iHopkins et"al]l2008l ). Thus, it is not surprising that we ob- 
serve such compact galaxies with high SFR. The models, 
indeed, predict a spheroidal- like morphology at the time 
of coalescence or just after, since the coalescence generally 
completes at about the sam e time as the gas first reaches 
the centre (|Cox et a l." 2008a). 

Moreover, our findings that a fraction of the SMGs have 
a compact morphology agrees with measurem ents of the ^as 
distr i bution from CO maps for many galaxies (^ Tacconi et al.1 
l200d l2008h . Therefore, these compact, highly star-forming 
systems are likely to be in the final phase of the merger, 
and are the transition link between starbursts and compact 
galaxies. 

A peculiar case is represented by the system 
SMMJ123635. 59+621424.1. We classify it as a merger event, 
as it appears as a double-nucleus system. However, its effec- 
tive radius and Sersic index are characteristic of those for 
compact galaxies. It is likely that we are measuring the com- 
pactness of the bulge component, given the faintness of the 
outer light. This systems also shows faint signs of potential 
spiral structure over a scale of 5 kpc. The SSFR of this sys- 
tem is the highest in our sample, compatible with being near 
the peak in star formation during the coalescence phase. It is 
worth to note that this object has the highest infrared lumi- 
nosity L = 10^^ °^ L© (from Michalowski, Hjorth, & W at^onl 
l2009h of our sample, hence it can be considered a HLIRG, 
supporting the picture where HLIRGs are galaxies in their 
maximal star formation periods triggered by interactions. 

In order to build an "illustrative" evolutionary se- 
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quence, we can order galaxies inside a given class accord- 
ing to the values of their SSFR. In Figure [H we show as a 
dashed line, the ev olution of the SSFR in a simulated galaxy 
merger, taken from iNaravanan et al ](l20ld). The simulation 
illustrates the evolution of the SFR in a major merger (with 
mass ratio 1:1) for a ^ 2 x 10^^ Mq dark matter halo. We 
have computed the SSFR taking the SFR of their Figure 
1 and dividing it by the final stellar mass (8 x 10^^ M©, 
which is roughly the maximum stellar masses of our sam- 
ple). The trend in SSFR shows a modest star formation 
rate in the pre-merger phase, in agreement with the value 
of our diffuse galaxies, followed by a steep rise during the 
starburst /merger, a peak in the coalescence phase, and then 
a rapid decline. The trend depicted by the simulations is 
matched by our data if we assume that the SSFR increases 
in the merging phase, and then declines for the compact 
galaxies. The peak is reached during the coalescence, as it is 
illustrated by the source SMMJ123635. 59+621424.1. In sim- 
ulations the diffuse/isolated systems are not fuelled by a new 
gas reservoir, as in the case for merging systems, and they 
simply exhaust the cold gas available. Hence their SSFR de- 
clines with time (see Cox et al. 2008b). Therefore, in the 
same plot, we have indicated by the coloured points the 
SSFR of our sample of SMGs, with the number in the x- 
axis indicating the position in the evolutionary sequence (in 
Figure [2] galaxies are ordered according to this sequence) . 
Note tha t we adopt the SFR derived from the IR luminos- 
ity from iMichalowski. Hiorth. WatsonI (| 20091 ). using the 
iKennicuttI (|l998l ) relation to computed the star formation 
rate. 
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The fact that we can put our galaxy sample in a SSFR- 
sequence that well matches that found in hydrodynamical 
simulations, strengthens the conclusion that the (morpho- 
logical) evolutionary sequence described above, where large 
diffuse systems transform themselves in compact remnants 
passing through a major merger event, is likely the formation 
mechanism for these galaxies into compact systems seen at 
slightly lower redshifts. Another piece of evidence support- 
ing this theoretical merging scenario is the size of the pro- 
genitor discs we have found. Theoretically t hese discs are 
expe c ted to have effective radius of 2-3 kpc (jPekel Coxl 
I2OO6I : iHopkins et al.""2009cl. This is in fact what our obser- 
vations show. Given the lack of NICMOS imaging for most 
of the sources and the insufficient resolution in the near-IR, 
we were forced to use ACS imaging to constrain our morpho- 
logical sequence. To better test such picture high-sensitivity 
imaging in near-IR would be required, and this could be 
achieved only when WFC3 imaging will be available. 



Summarising, our data can not reject the evolutionary 
picture depicted by theoretical models in which the precur- 
sors of the superdense galaxies are massive, gas-rich discs 
at z ^ 2 — 3, which evolve into compact remnants through 
dissipative major mergers. Moreover, a cold gas accretion 
driven scenario for the formati on of the compact massive 
galaxies, as the one proposed bv lPekel et al.l (|2009l l perhaps 
can not be easily supported by our data, since in this case 
we would likely not observe such a large diversity of size and 
structure for these progenitors as we observe here. 



Table 1. Properties of SMGs 



ACS NICMOS 



Galaxy^ 






log- _M_ c 


SSFR_£ 

[Gyr-i] 


Rab^ 


Re 

arcsec 


Re 

[kpc] 


n 


b/a 


Re 

arcsec 


Re 

[kpc] 




n 


b/a 


Class £ 


SMMJ123711 


98+621325.7 


1.992 


10.99 


1.72 


25.8 


0.36 


3.03+0.03 


0.69+0.03 


0.79+0.00 £ 












D 


SMMJ123606 


85+621021.4 


2.509 


11.69 


1.17 


25.3 


0.25 


2.03+0.04 


2.38+0.00 


0.70+0.00 












D 


SMMJ123618 


33+621550.5 


1.865 


11.40 


0.66 


25.9 


0.34 


2.84+0.12 


2.19+0.14 


0.66+0.12 












D 


SMMJ123707 


21+621408.1 


2.484 


11.65 


0.38 


26.0 


0.43 


3.46+0.39 


2.54+0.73 


0.43+0.11 












D 


SMMJ123712 


05+621212.3 


2.914 


11.80 


0.20 


25.5 


0.25 


1.97+0.04 


0.98+0.02 


0.67+0.00 


0.21 


1.60+0 


11 


1.51+0.36 


0.72+0.09 


D 


SMMJ123616 


15+621513.7 


2.578 


11.67 


0.80 


25.7 


0.10 


0.83+0.05 


2.64+0.76 


0.74+0.07 


0.02 


0.16+0 


12 


3.91+2.69 


0.47+0.43 


M 














0.07 


0.53+0.04 


0.85+0.30 


0.53+0.05 


























0.11 


0.86+0.07 


0.69+0.13 


0.66+0.06 


0.28 


2.23+0 


17 


1.00+0.00 


0.97+0.05 




SMMJ123622 


65+621629.7 


2.466 


11.62 


1.72 


25.4 


0.09 


0.74+0.04 


1.50+0.02 


0.12+0.01 


0.35 


2.84+0 


07 


0.37+0.02 


0.41+0.03 


M 














0.58 


4.66+0.43 


1.96+0.17 


0.10+0.00 


0.54 


4.38+0 


73 


2.23+0.04 


0.35+0.04 




SMMJ123553 


26+621337.7 


2.098 


11.03 


5.73 


24.7 


0.43 


3.61+0.15 


2.75+0.12 


0.28+0.00 












M 














0.14 


1.12+0.06 


0.51+0.04 


0.22+0.00 














SMMJ123635 


59+621424.1 


2.005 


11.08 


7.23 


24.2 


0.17 


1.39+0.28 


0.87+0.00 


0.78+0.00 


0.44 


3.64+0 


12 


1.97+0.08 


0.81+0.04 


M 


SMMJ123632 


61+620800.1 


1.993 


11.03 


4.52 


23.6 


0.01 


0.08+0.04 


4.65+1.22 


1.00+0.00 


0.05 


0.39+0 


73 


2.78+3.33 


0.45+0.42 


C 


SMMJ123600 


15+621047.2 


1.994 


11.34 


2.53 


25.1 


0.06 


0.53+0.10 


3.57+0.66 


1.00+0.00 


0.19 


1.60+0 


18 


0.63+0.28 


1.00+0.00 


C 


SMMJ123606 


72+621550.7 


2.416 


11.15 


1.27 


23.6 


0.002 


0.02+0.02 


4.47+3.12 


1.00+0.00 












c 



^ Galaxy ID, objects without ID are the subcomponents in the merging systems. 
^ Spectroscopic red shift fromlChapman et al.l (120051) . 

^ Stellar mass from lMichalowski. Hiorth. Sz WatsonI (|2009l ) converted to a Kroupa IMF. 
^ Specific star formation rate from lMichalowski. Hiorth. &: WatsonI (l2009l ). 
^ Magnitudes in the R-band (AB system") fromlC hapman et al.l (|2005f i 

f Morphological class: D stays for disc-like objects, M for ongoing mergers and C for compact galaxies (see Sect. 4). 
9 Measurements with errors equal to 0.00 mean that the parameters have been fixed in the fitting procedure. 
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